Reversible Formation of a Bose-Einstein Condensate 
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We present a method of adiabatically changing the local phase-space density of an ultracold gas, 
using a combination of magnetic and optical forces. Applying this method, we observe phase-space 
density increases in a gas of sodium atoms by as much as 50-fold. The transition to Bose-Einstein 
condensation was crossed reversibly, attaining condensate fractions of up to 30%. Measurements of 
the condensate fraction reveal its reduction due to interactions. 



The physical properties of atomic gases change dra- 
matically when quantum degeneracy is reached, i.e. when 
the ground state population approaches unity Q. Re- 
cent successes in reaching quantum degeneracy with 
Bose gases [|f-fi} have relied on non-adiabatic, irreversible 
methods such as laser and evaporative cooling. The pos- 
sibility of changing the ground state population by an 
adiabatic change in the trapping potential had been over- 
looked for quite some time ||. Indeed, in the case of an 
ideal gas, adiabatic changes in the strength of the trap- 
ping potential do not change the ground state popula- 
tion HQ ■ However, Pinkse and collaborators || recently 
showed, both theoretically and experimentally, that by 
changing the form of the trapping potential, the popu- 
lation in the ground state can be changed adiabatically. 
For a non-degenerate gas the ground state population is 
identical to the phase-space density T = nA T , where n is 
the density of the gas and At is the thermal de Broglie 
wavelength. Within the type of trap deformations con- 
sidered in Rcf. [|| the maximum increase of phase-space 
density is limited to a factor of 20. 

In this Letter, we show that a more general deforma- 
tion of the trapping potential can increase the phase- 
space density by an arbitrary factor, and we describe 
an implementation of this scheme using a combination 
of magnetic and optical forces. Furthermore, we demon- 
strate the ability to cross the Bose-Einstein condensation 
(BEC) phase transition and to change the condensate 
fraction of a Bose gas in a reversible way. 

Adiabatic increase in phase-space density. The 
type of trap deformations which we study can be un- 
derstood with the following "two-box" model. Consider 
a classical gas of N atoms confined in a box of volume 
Vq = V1 + V2 with an initial phase-space density Tq. Sup- 
pose that the potential within a sub-volume V2 of the 
box is lowered to a final well-depth U. In this final po- 
tential, the gas equilibrates at a temperature T, and the 
density in V% will be higher than that in Vi by the Boltz- 
mann factor e c/ /' CBT . Using the condition of adiabaticity 
and constant particle number, one obtains the relative 
increase of phase-space density in V2 compared to that 
in Vq before compression: 



in(r 2 /r ) = 
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For deep potential wells, where U/ksT » ln(Fi/p2), 
there is no increase in phase-space density since all of the 
gas becomes confined in V2, and the adiabatic deforma- 
tion corresponds simply to a uniform compression of the 
gas. For shallow potential wells (U/k B T < ln^i/V^)), 
the phase-space density in V2 increases as e u / kBT . As 
U is varied between these limits, the phase-space den- 
sity increase reaches a maximum which is greater than 
(V1/V2) 1 / 2 . Thus, by choosing an extreme ratio of vol- 
umes V1/V2, an arbitrarily large increase in phase-space 
density r2/ro is possible. 

To demonstrate this phase-space density increase in 
a gas of trapped atoms, a narrow potential well (anal- 
ogous to V2) was added to a broad harmonic potential 
(corresponding to Vi ) by focusing a single infrared laser- 
beam at the center of a magnetic trap. First a gas of 
atomic sodium was evaporatively cooled to a tempera- 
ture higher than the BEC phase transition temperture 
in the cloverleaf magnetic trap |J . The number of atoms 
and their temperature were adjusted by varying the final 
radio frequency (rf) used in the rf-evaporative cooling 
stage JlOt] . Afterwards, the magnetically trapped cloud 
was decompressed by slowly reducing the currents in the 
magnetic trapping coils. Time-of-flight absorption imag- 
ing was used to characterize the cloud. The total num- 
ber of atoms N was determined by integrating the col- 
umn density across the cloud, and T was determined by 
one-dimensional Gaussian fits to the wings of the den- 
sity distribution. From these, we determined the fugac- 
ity z of the gas by the relation gz(z) = N(hui/kBT) 3 , 

and then its phase-space density by r = 33/2(2), where 
00 

g n [z) — z% /i n PJl- Here, Q is the geometric-mean 

i=l 

trapping frequency of the magnetic trap, as determined 
by in situ measurements [|12| . These phase-space density 
measurements were calibrated with images from magnet- 
ically trapped clouds at the phase transition. 

The optical setup was similar to that used in Ref . jl3) . 
The infrared laser power was gradually ramped-up from 
zero to a power P c at which the onset of BEC was seen 
in time-of-flight images of clouds released from the de- 
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formed trap; this implied Tf = 33/2(1) = 2.612 for the 
final phase-space density. The depth of the optical poten- 
tial well was given by U c /kB = 37 [iK P c /wq (^m 2 /mW), 
where wo is the 1/e 2 beam- waist radius at the focus. The 
ramp-up time was made long enough to ensure that the 
trap deformation was adiabatic, but also short enough to 
minimize heating and trap loss. Ramp- up times of up to 
10 s were used. 
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FIG. 1. Phase-space density increase to reach BEC vs. 
normalized well-depth (a). Various trap settings were used: 
uj = 2tv x 100 Hz, wo — 9/im (triangles); uj — 2ty x 100 Hz, 
Wo — 18 fim (open circles) ; and uj — 2ir x 33 Hz, wo = 18 /im 
(closed circles). The solid line gives the prediction of Eq. (Q). 
(b) Prior to deformation, the harmonic trapping potential 
V{r) holds a cloud above BEC the transition temperature 
(jj, < 0). (c) When a potential well with depth U = —fx is 
added, a small condensate forms. 

Fig. |l|a shows the increases in maximum phase-space 
density which were measured at three different settings 
of the trap parameters. A maximum increase by a factor 
of 50 was obtained. Condensates were observed in clouds 
with temperatures as high as 5 fiK. Further increases 
were hindered by limitations in laser power and by limits 
to the ramp-up time set by the various heating and loss 
processes in the deformed trap. The scatter in the data 
is primarily due to statistical errors in our measurements 
of Uc/ksT, at the level of 30%, which arise from mea- 
surements of P c and wo , and in the determination of the 
transition point. 

The well-depth U c required to reach BEC can be un- 
derstood by a simple model depicted in Fig. |l]. We begin 
with a gas in a harmonic trap above the BEC transition 
temperature, i.e. its chemical potential \i < (Fig. |]b). 
Its initial phase-space density is given by r = 53/2(2), 
where z — e t 1 / k B r r j s ^he f U g ac ity. In Fig. Qc, a narrow 
potential with depth U c is added until U c = —fi, at which 
point the gas begins to Bose condense. The phase-space 
density thus reaches the critical value 33/2(1), so that the 
increase in maximum phase-space density is given by 

Tf = ff3/ 2 (l) (2) 

r o g 3/2 (exp(-U c /k B T)) 

This prediction, shown in Fig. [l]a, describes our data well, 
and accounts for the universal behaviour of our measure- 
ments over a wide range of temperatures and well-depths. 
Note that in this simple model we made the implicit 



assumption that the initial and final temperatures of the 
cloud were equal. We can remove this assumption by 
considering instead the condition of constant entropy. 
The entropy per non-condensed particle in a harmon- 
ically confined Bose gas is determined uniquely by its 
fugacity z @: 

S , , , qa(z) 
( z )=4^-]nz . 3 

This equation describes the entropy of the gas before 
compression, with the fugacity given by zq. After com- 
pression, because of the small volume of the poten- 
tial well, the entropy per particle is approximately that 
of a harmonically trapped gas (Eq. (||)) with fugacity 
Zf = e~ u "/ kBT . Here T is the final temperature of the 
gas, which is generally higher than the initial tempera- 
ture jL5|. Constant entropy then implies zq — Zf, and 
thus one obtains Eq. (||) where T is now the final tem- 
perature, as we have plotted in Fig. |l|a. 

One may also consider the process of adiabatically in- 
creasing the phase-space density as a change in the den- 
sity of states D(e) of the system. By increasing the well- 
depth in a small region of the trap, we lowered only the 
energy of the ground state and a few excited states. Thus 
r, a local quantity, is maximally increased as the ground 
state energy is brought ever closer to the chemical po- 
tential, while the entropy, a global property of the gas, is 
unchanged by the minimal modification of D(e). 
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FIG. 2. Momentum distributions of the thermal cloud 
at the BEC phase transition show no Bose-enhancement for 
the deformed trap (a), but a clear Bose-enhancement for the 
purely magnetic trap (b). Both distributions show a small 
condensate peak. Lines show a Maxwell-Boltzmann distri- 
bution (dotted) and a Bose-Einstein distribution (for z = 1) 
(dashed) for clouds in a harmonic trap, which were fit to the 
thermal wings. The momentum distribution is shown as a 
profile across absorption images taken after 40 ms of ballistic 
expansion. 

The fact that global properties of the gas are not af- 
fected by the trap deformation can also be seen in the 
momentum distributions probed by timc-of-flight imag- 
ing [g|Jl^]. The onset of BEC in the combined opti- 
cal and magnetic trap is signaled only by the formation 
of a condensate peak. The remaining thermal cloud is 
well fit by a Maxwell-Boltzmann distribution, which de- 
scribes a magnetically trapped cloud far from condensa- 
tion (Fig. 0a). In contrast, at the BEC transition in the 
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harmonic magnetic trap, the momentum distribution of 
the thermal cloud is clearly Bose-enhanced at low mo- 
menta (Fig. |b). 

Adiabatic condensation. We now turn to the stud- 
ies of adiabatic, i.e., reversible, condensate formation. A 
cloud of about 50 x 10 6 atoms was evaporatively cooled 
down to the the transition temperature in the mag- 
netic trap, at trap frequencies of cu r = 2n x 20 Hz and 
uj z — 2tt x 13 Hz in the radial and axial direction, respec- 
tively. The power of the infrared laser beam (of radius 
wq = 20 jtzm) , was ramped up over 1 s and held at a 
constant power for a dwell time of 1.5 s. Condensate 
fractions as small as 1% could be distinguished from the 
normal fraction by their anisotropic expansion in time- 
of- flight images [0,0. The condensate number No was 
determined by subtracting out the thermal cloud back- 
ground using Gaussian fits to the thermal cloud in regions 
where the condensate was clearly absent. 
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FIG. 3. Condensate formation by adiabatic trap deforma- 
tion. The condensate fraction is plotted against the (nor- 
malized) well-depth. Lines show predictions for an ideal gas 
(solid) and for an interacting gas (dashed). 

As shown in Fig. 3, the adiabatic trap deformation 
yielded condensate fractions of up to 15%; we observed 
condensate fractions of 30% with different trap settings. 
By varying the dwell time, we confirmed that clouds 
for which U /ksT < 1.5 suffered no significant losses of 
condensate number due to three-body decay or heating, 
whereas those points with higher values of U /ksT were 
affected by such losses. 

For an ideal Bose gas, the result of this adiabatic 
change can be understood as follows. Before compres- 
sion, the cloud of N particles at the BEC transition has 
an entropy Si given by Eq. (||) as Si — N x 45 4 (l)/g 3 (l). 
After compression, the situation is similar to that indi- 
cated in Fig. lc, i.e. because of the small volume of the 
attractive well, the cloud is well-described as a harmoni- 
cally trapped gas with /i = —U. Thus, accounting for the 
fact that condensate particles carry no entropy and again 
using Eq. (^|), which gives the entropy per non-condensed 
particle, we equate the entropy before and after compres- 
sion and obtain 

= 4 g4 (l)/ g3 (l) 

N A gi {e- u / k BT)/ g ^ e -u/k B T) + u / kBT - \ > 



However, this simple prediction does not describe our 
findings well. The theory described above has two short- 
comings. First, the approximation of using Eq. (^) for 
the deformed trap is not strictly valid. However, cal- 
culations which accounted for the true shape of the de- 
formed potential changed the prediction of Eq. (Q) only 
for U/ksT > 1, and only slightly improved the fit to our 
data. 

A second shortcoming is the neglect of interactions. It 
has been shown that in harmonic potentials, the conden- 
sate fraction in a Bose gas with repulsive interactions is 
reduced in comparison to that of an ideal gas |lj],|l8| ■ To 
estimate this effect in our system, we use the "semi-ideal" 
model of Ref. |l8) . The thermal cloud is described as an 
ideal gas for which the chemical potential is raised by 
gno = 4Trh 2 ano/m, where no is the maximum conden- 
sate density, m the mass of sodium, and a — 2.75 nm its 
scattering length |jl9). This simply corresponds to using 
Eq. (jij) with the substitution U — » U — giiQ. We de- 
termined no using the Gross-Pitaevskii equation in the 
Thomas- Fermi limit [pof , and a harmonic approximation 
for the deformed trap potential at its center. 

This approach predicts a significant reduction of the 
condensate fraction (Fig. 3, dashed line), and the im- 
proved agreement with our data is strong evidence for 
this effect. In contrast to related studies in purely har- 
monic traps [pyi6|], which did not show evidence for in- 
teraction effects, this depletion is strongly enhanced by 
the shape of the potential we are using. The mean-field 
energy of the condensate gno is large because the con- 
densate forms in the tight optical potential, while the 
transition temperature T c is small since it is determined 
by the weak magnetic potential. 

The reversibility of crossing the BEC phase transition 
was demonstrated by preparing a magnetically trapped 
cloud just above T c . We then sinusoidally modulated 
the power of the infrared light at 1 Hz, between zero 
and 7 mW. This modulation frequency was significantly 
smaller than the magnetic trap frequencies (uj r — 2tt x 
48 Hz and lj z = 2ir x 16 Hz). These low frequencies and a 
large optical focus (wo — 18 fim) were used to minimize 
trap loss due to inelastic collisions. 

During the first seven condensation cycles the conden- 
sate fraction oscillated between zero and 6% (Fig. 4); 
later probing showed repeated condensation for at least 
15 cycles. The peak of these oscillations decreased slowly 
in time. The temperature also oscillated, with an ampli- 
tude of about 100 nK, while gradually rising by about 
10 nK/s. This heating and the decrease of the peak con- 
densate fraction are consistent with similar behaviour in 
clouds held at a constant infrared power, which result 
from beam jitter, spontaneous scattering, and three-body 
decay Jl3| . Thus, within the stability limitations of our 
optical setup, the repeated crossing of the BEC phase 
transition appears fully adiabatic. 
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FIG. 4. Adiabatic cycling through the phase transition. 
Shown is the condensate fraction (a) and the temperature 
(b) vs time for the case of a modulated infrared beam (closed 
circles), an infrared beam ramped up to a constant power 
(open triangles), and no infrared light (open squares). The 
solid lines are guides to the eye. 

This method of creating condensates provides insight 
into their formation, which was recently studied experi- 
mentally and theoretically |22| . For example, in the 
experiment described above, the condensate fraction was 
found to lag about 70 ms behind the modulation of the 
laser power, which is a measure for the formation time. 

In other experiments, by switching on the infrared light 
instantly, we observed condensation on timescales much 
faster than the oscillation periods in the magnetic trap 
and along the weakly confining axis of the optical trap. 
The resulting condensates showed striations in time-of- 
night images, indicating that the condensates formed into 
excited states of the deformed potential. Such studies of 
shock-condensation might give new insight into the for- 
mation of quasi-condensates and condensation into ex- 
cited states p3p^l. 

In conclusion we have demonstrated the adiabatic 
Bose-Einstein condensation of an ultracold gas of atomic 
sodium. Changes in the trapping potential resulted in 
large phase-space density increases and allowed for re- 
peated crossings of the BEC phase transition. This 
method allows for detailed studies of condensate forma- 
tion and the phase transition. The combined trapping 
potential widens the range of trap parameters over which 
BEC can be studied. This was used to strikingly enhance 
the role of interactions, and led to higher transition tem- 
peratures (up to 5 /iK) than achieved in purely magnetic 
traps. 

We thank Michael R. Andrews for discussions. This 
work was supported by the Office of Naval Re- 
search, NSF, Joint Services Electronics Program (ARO), 
NASA, and the David and Lucile Packard Foundation. 
A. P. C. acknowledges additional support from the NSF, 
D. M. S. -K. from JSEP, and J. S. from the Alexander 
von Humboldt-Foundation. 



A. Griffin, D.W. Snoke, and S. Stringari, Bose-Einstein 
Condensation (Cambridge University Press, Cambridge, 
1995); Georgia Southern Universit y BEC home page, 



[2] 

[3] 



[5] 
[6] 
[7] 



[8] 
[9] 
[10] 



[11] 

[12] 
[13] 
[14] 
[15] 



[16] 
[17] 



[18] 

[19] 

[20] 
[21] 
[22] 

[23] 

[24] 



http://amo.phy.gasou.edu/bec.html 



M.H. Anderson et al, Science 269, 198 (1995). 
K.B. Davis et al, Phys. Rev. Lett. 75, 3969 (1995). 
C.C Bradley, C.A. Sackett, and R.G. Hulet, Phys. Rev. 
Lett. 78, 985 (1997), see also: C.C.Bradley et al, Phys. 
Rev. Lett. 75, 1687 (1995). 

W. Ketterle and D.E. Pritchard, Phys. Rev. A 46, 4051 
(1992). 

F. Reif, Fundamentals of Statistical and Thermal Physics 

(McGraw-Hill Book Company, New York, 1965) 

For weakly interacting Bose gases this argument has to 

be modified; see: M. Houbiers, H.T.C Stoof, and E.A. 

Cornell, Phys. Rev. A 56, 2041 (1997). 

P. W. H. Pinkse et al, Phys. Rev. Lett. 78, 990 (1997). 

M.-O. Mewes et al, Phys. Rev. Lett. 77, 416 (1996). 

W. Ketterle and N.J. van Druten, in Advances in Atomic, 

Molecular, and Optical Physics, edited by B. Bederson 

and H. Walther (Academic Press, San Diego, 1996), Vol. 

37, p 181. 

This assumes an ideal Bose gas and harmonic confine- 
ment; see: V. Bagnato, D. E. Pritchard, and D. Klepp- 
ner, Phys. Rev. A 35, 4354 (1987). 



D. M. Stamper-Kurn et al., preprint cond-mat/9801262 
(unpublished). 

D. M. Stamper-Kurn et al., Phys. Rev. Lett. 80, 2027 
(1998). 

K. Huang, Statistical Mechanics, second edition (Wiley, 
New York, 1987). 

This approximation is valid for shallow potential wells 
(U/ksT <C ln(Vi/V2)), as discussed in the two-box 
model. Numerical calculations confirmed its validity for 
all data in Fig. [j]a. Note that the effective volume of the 
magnetic trap (Vi) increases with temperature. 
J. R. Ensher et al, Phys. Rev. Lett. 77, 4984 (1996). 
V.V. Goldman, I.F. Silvera, and A.J. Leggett, Phys. Rev. 
B 24,2870 (1981); R.J. Dodd, K. Burnett, M. Edwards, 
and C.W. Clark, Acta Physica Polonica 93,45 (1998); S. 
Giorgini, LP. Pitaevskii and S. Stringari, Phys. Rev. A 
54, R4633 (1996). 

M. Naraschewski and D.M. Stamper-Kurn, preprint 



cond-mat/9803216 (unpublished). 

E. Tiesinga et al, J. Res. Natl. Inst. Stand. Technol. 101, 
505 (1996). 

G. Baym and C.J. Pethick, Phys. Rev. Lett 76, 6 (1996). 

H. -J. Miesner et al, Science 279, 1005 (1998). 

C.W. Gardiner, P. Zoller, R.J. Ballagh, and M.J. Davis, 
Phys. Rev. Lett. 79, 1793 (1997). 

Yu. Kagan, B.V. Svistunov, and G.V. Shlyapnkov, Sov. 
Phys. JETP 75, 387 (1992). 



C.W. Gardiner et al, preprint cond-mat/9801027 (un- 
published) . 



4 



